Surface States of Topological Insulators 
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We develop an effective bulk model with a topological boundary condition to study the surface 
states of topological insulators. We find that the Dirac point energy, the band curvature and the 
spin texture of surface states are crystal face-dependent. For a given face on a sphere, the Dirac 
point energy is determined by the bulk physics that breaks p-h symmetry in the surface normal 
direction and is tunable by surface potentials that preserve T symmetry. Constant energy contours 
near the Dirac point are ellipses with spin textures that are helical on the S/N pole, collapsed to 
one dimension on any side face, and tilted out-of-plane otherwise. Our findings identify a route to 
engineering the Dirac point physics on the surfaces of real materials. 

PACS numbers: 71.70.Ej, 73.20.-r, 73.22.Dj 



Introduction. — The discovery [Ij-iG, of topological insu- 
lators (TI's) and the synthesis [7HTU] of three dimensional 
materials that realize their physics has opened up a new 
field in solid state physics. Particular interest has focused 
on the TI surface states with point degeneracies that 
are topologically protected by time reversal symmetry 
(T) when a trivial insulator with topological index[Tl [2] 
Z2 = 1 is joined to a TI with Z2 = —1. On the cleav- 
age surface of Bi2Se3, a TI with a single Dirac cone 
has been identified [,.9 f,l_5j by angle resolved photoemission 
(ARPES) and scanning tunneling microscope (STM) ex- 
periments. The surface states form a spin-momentum 
locked helical metal with conduction and valence bands 
exhibiting opposite helicities near the Dirac point, and 
develop hexagonal warping [101 [TSHTT] on top of the cir- 
cular constant energy contour away from the Dirac point. 

Interest in the topological surface bands has led to 
the development of various bulk continuum models with 
the surface termination described by ad-hoc fixed node 
[H [T7H^ or by "natural" boundary conditions pi] . 
where the existence of a topological surface state and its 
Dirac point energy are treated as inputs to the theories. 
The choice of boundary condition has profound conse- 
quences for the spatial form of the surface state wave- 
functions and their interactions with external fields and 
absorbed species. In this work we derive the appropriate 
topological boundary condition (TBC) at the surface from 
the matching of bulk waves to evanescent vacuum states. 
As expected a topologically protected Dirac point occurs 
in this model but its energy position and band curvature 
are both crystal face-dependent owing to bulk terms that 
break particle-hole (p-h) symmetry. Furthermore, the in- 
teraction with surface-localized potentials that preserve 
T symmetry shifts the energy position of the Dirac de- 
generacy and provides a robust degree of freedom for tun- 
ing the surface state spectrum into a convenient energy 
range to exploit their topological properties. The associ- 
ated spin texture is determined by the bulk symmetries 
and depends on the crystal face angle. Unlike the helical 
metal on a cleavage surface, the spin texture near Dirac 
point is compressed to a single direction on a side face, 



and is tilted out-of-plane otherwise. 

We start from a description of the low-energy model 
of Bi2Se3, which applies generally to other TI's with the 
same crystal structure with space group R3fh. At the 
origin of Bi2Se3 Brillouin zone F, the effective Hilbert 
space near the bulk gap is spanned by states with angular 
momentum mj = zt^ and parity V — H. Because of the 
spin-orbit coupling (SOC), \pz t) and |pz |) states are 
mixed with |p+ |) and |p_ f) states, respectively. Since 
the crystal-field splitting is much larger than the SOC, Pz 
orbitals dominate and mj pseudospin is proportional to 
the electron spin. The P = ±1 hybridized states can be 
labelled approximately by |-|-) state from Bi atoms while 
|— ) from Se, due to the large energy difference between 
Ap (Se) and 6p (Bi) principal quantum levels. 

Besides T and the parity inversion (V) symmetries, 
Bi2Se3 crystal structure has threefold rotational (C3) 
symmetry along the z perpendicular to the quintuple 
layers, and twofold rotational (C2) symmetry along TM 
direction. By convention we choose the parity opera- 
tor V — T2 and the time reversal operator T — iKuy 
where K is the complex conjugate operation. Therefore, 
to quadratic order in k the k ■ p bulk Hamiltonian that 
preserves the above four symmetries has a unique form 

H = {co + c^kl + c^^kl) + {—nif) + m^kl + m„kl)Tz 

+ WzfczTy -|-W,|(fcyCrx- fcxO'y)Tx, (1) 

where we assume mz,m„,Vz,v„ > and h= 1 hereafter. 
The first parentheses is a scalar term that preserves T 
and V symmetries but breaks the p-h symmetry. The 
quadratic scalar terms intrinsically give rise to the cur- 
vature of Dirac surface bands and the rigid shift of Dirac 
point from the middle of the bulk gap, as we will demon- 
strate in the following. Eq.(fl]) has both T and V sym- 
metries with a topological index Z2 = sgn(— toq). 

Topological boundary condition. — To investigate the 
surface states, we use a TBC in which the mass term 
changes sign across the surface. On the TI side, toq = m 
is positive and half of the bulk gap at fc = 0, and cq — 
to define the middle of the A; = gap as energy zero. On 
the vacuum side, i.e., a trivial insulator with an infinite 
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FIG. 1: (Color online) The definition of an arbitrary face of 
a Bi2Se3-type topological insulator. The south (S) pole de- 
notes the (001) surface parallel to the quintuple layers and the 
equator represents the side faces perpendicular to the quintu- 
ple layers. For any face E(6l), fcs _L E and ^2 ~ ky. 



gap mo = —M where M -^ -f-cx), both c^ and c„ vanish 
since the vacuum has p-h symmetry around cq. 

We first focus on the one dimensional quantum me- 
chanics of (001) surface turning off the quadratic terms 
in Eq.([I]). Notice that the two spin flavors are decoupled 
for this eigenvalue problem. The spin f states in {r^, a} 
representation are determined by the following coupled 
Dirac equations 






m[z) — E 



V'2 



0. 



(2) 



The wave function ■0^ = {'ipi,'ip2y is continuous across 
the surface, integrating Eq.(l2]) over the vicinity of the 
surface. This continuity condition leads to a nontrivial 
solution which is isolated in the middle of the bulk gap 
and is localized on the surface. This midgap {E — 0) 
state has a simple and elegant exact solution, 



'ipk„.,tix,y,z) 



J{ky,x+kyy) 



e"-'^^- • ■^^"'0{z) 



e-«^ z > (TI) 
e'"«^ z < (Vac) ' 



(3) 



where k — m/vz, kq — M/v^ and A is a normalization 
factor. (/)(z) is evanescent on both sides of the surface 
where the mass mo changes sign, analogous to Jackiw 
and Rebbi solution[55] of a two-band Dirac model. The 
spin l solution can be obtained by ipi = Ttp^- This 
isolated midgap state at fci, = is identified as the Dirac 
point and at finite fc,, spreads into a perfect Dirac cone- 
the ideal topologically protected surface bands. 

Surface state spin texture. — On the (001) surface, the 
midgap solution of the boundary problem at fc,, = is 
determined by the operators t and is free under any ro- 
tation of the operators cr. This t ^ cr algebra of the 



boundary problem can be generalized to any surface with 
r replaced by Si and cr by 52- For an arbitrary crystal 
face E(6') defined in Fig{l] the algebraic structure is 



Si 

S2 






- /3cryTy, aTy - /JO-yTx, T^} , 

■ /3crzTz, o-y, acTz + /3a-xTz} ., 



(4) 

Uz COS O/v^ 

These new pseudospins satisfy 

^in^h\ ~ 2i<5abe''''^'S'a. Rewritten in this new pseudospin 



where V3 = \/{vz cos 9)^ + (wi, sin 9)"^ 
and /3 = w,, sin 0/^3 

basis, Eq.([T]) reads 



n 



- ruQ SI + (v^ks + voh) Sf 
+ {v,ky S^ - viki Si) Sf , 



(5) 



where Vq = (v'^^ — v^)sm9 cos9/v3 and vi = v^v^/v^. By 
matching the eigensystems of TI and vacuum sides, we 
obtain the S(0) surface states similar to Eq.(|3| where 
T and <T are replaced by Si and ^2, respectively. Note 
that K = m/wa -t- ivoki/v^ in general, fei coupling to 
S\ neither influences the energy spectrum nor the spin 
texture, only making the evanescent states oscillate and 
giving a phase accumulation along ki away from Dirac 
point. Therefore, ignoring the quadratic corrections, we 
first derive the effective surface state Hamiltonian for an 
arbitrary face S(0) to the linear order. 



n'-^\9) = V,ky S^ 



viki Si , 



(6) 



from which we can further explicitly demonstrate the sur- 
face state spin texture on T,{9): 

I , , V2V„kyC0s9 



^y)e 



± 



vzJvlkl ^vlk^ 

-VzVyki 
vsJvjkl +vf,k^ 



0, 



(7) 



where + (— ) denotes the conduction (valence) band. 
The electron real spin (s) is proportional to but always 
smaller than {cr) due to the SOCJ23]. In the local coordi- 
nates, (tTi)e = {<7x)ecos9 and {a3)g = ((Jx)esin0. Eq.([7]) 
indicates that the surface state spin texture is rather dif- 
ferent from face to face while its pseudospin (S'2) has a 
universal structure on the elliptic constant energy con- 
tour near the Dirac point. Clearly, there is no spherical 
symmetry to guarantee that the spin and orbital struc- 
tures are the same anywhere on a TI sphere. The crystal 
face-dependent spin texture is helical on the south and 
north poles, is compressed to a single dimension along 
the equator, and is tilted out-of-plane otherwise. The 
surface state anisotropy and spin textures for different 
faces are compared in FigJ2J 

The surface state spin textures on the poles and the 
equator of a TI sphere can be understood by symme- 
tries. For (100) face, the C2 symmetry along TM (x) di- 
rection forbids Ux coupling to any in-plane momentum. 



In linear order C3 symmetry upgrades to continuous ro- 
tational symmetry, consequently, the surface normal spin 
is zero along the TI equator. On the two poles, the mir- 
ror symmetry and the C3 symmetry insure that the spin- 
momentum locking into the form of kya^ — ky^cjy. a^ 
decouples to any momentum in the linear order, taking 
into account all the four symmetries. 

It's interesting to point out that the surface band is 
the positive eigenstate of S'f and the chiral counterpart 
is separated and localized on the opposite face. Thus 
the surface state Hilbert space is reduced by half and 
this pseudospin polarity (or chirality) blocks the back 
scattering on the surface as a result of interplay between 
T symmetry and band inversion physics with TBC. 

Dirac point energy and surface potentials. — The ab- 
sence of spherical symmetry in the bulk requires that 
the Dirac point has different energies on different crystal 
faces. When the quadratic mass terms and p-h symmetry 
breaking terms in Eq. M are turned on, the surface state 
wave function and spectrum are changed. In the coupled 
Schrodinger-type Eq.Q, the components of wave func- 
tion and their slopes are all continuous across the bound- 
ary. We find that the effective Hamiltonian for face S(0) 
have two important [23] corrections: 



n^^^e) = c,kl + {c.^sm'^9 + c„cos^0)kj, (8) 

^2 , „ „;„2, 



H^^(0) = 



y 

Cz cos" 



+ C|| sm 



cos^ 9 + m,, sin 



m. 



(9) 



The breaking p-h symmetry terms give the surface Dirac 
cone a parabolic curvature described by Eq. ^ and shifts 
the Dirac point from the midgap to a nonzero energy 
given by Eq.Q. These two effects help to explain the 
origin of nonzero Dirac point energies and nonlinear 
Dirac cones of cleavage surface states observed in ARPES 
experiments UHl [15H23 • 

Although the surface state solution obtained by Eq.([2]) 
with TBC remains topologically stable, it is essential 
to understand how localized surface potentials influence 
the surface spectra and the associated wave functions as 
well. We focus on potentials that preserve T symmetry. 
Among these six types: /, r^, Tx and aTy list in Tablelll r^ 
and (7Ty break V symmetry. It turns out that the same 
potential could play different roles on different crystal 
faces, as shown in Table |l] 

On the south pole of a TI sphere, a potential AqS^z) ■ 
2vz/m-I changes the wave function continuity conditions, 
integrating Eq.S including surface potentials over the 
vicinity of z = 0. The amplitudes of ■01 and ip2 are either 
weakened or enhanced across the surface with changes are 
always opposite to each other. Consequently, to match 
the evanescent solutions on the vacuum and the mat- 
ter sides, the surface band energies are rigidly shifted by 
(JiS'dp- Tyi potentials have similar effects to the /-type 
although Tx breaks V symmetry. For a r^ potential, it 
simply modifies the mass term on the surface and it is 
not surprising that it does not affect the surface spectra. 
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FIG. 2: (Color online) Dirac cones and spin textures of sur- 
face states on faces with (a) 9 — 0, (b) 6 — 7r/4 and (c) 
9 — n/2. The equal-energy contours are from —80 meV to 
160 meV with 10 meV increment relative to the Dirac point. 
Surface band curvatures are taken into account. All three 
panels are in the same scale and the parameters are adopted 
from DFT results. The lower panels only show the spin tex- 
tures at 160 meV in the fci — fcy planes where fcy is the vertical 



A Tz potential tunes the amplitudes of ■0i and ^2 in the 
same manner on the matter side but this gives no observ- 
able effect since "01 (-z) = 02 (^) and "01.2(2 < 0) ^ are 
still valid. Similarly, a potential like CTnTy (AuCTh = A • ct) 
does nothing to the surface spectra and 0i/02; however. 



TABLE I: Summary of the influence of T symmetry aUowed 
momentum-independent surface potentials A5{r3) ■ 2v2./m on 
the inversion symmetry, and the wave function continuity and 
the Dirac point energy of surface states. This A represents 
different surface potentials and their corresponding energy 
scales: Aq/, AmTz, AeS^ and AnS'2S'i. For the a^Ty col- 
umn, only the results for 55 = 1 state is shown and their 
complex conjugates represent the results for 5'2 = — 1 state. 
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it couples the two spin flavors and shifts their phase in 
an opposite way. Note that 4'i = V'2 is always true on 
the vacuum side since M —)' oo. The above results are 
summarized in Table IT] providing sufficient information 
to construct the surface state wave function and to engi- 
neer the Dirac point energy position. 

On an arbitrary face S(^), the types of surface poten- 
tials are the same but their combinations and the corre- 
sponding roles are rearranged. This can be fully under- 
stood by the fact that the spin-orbital structure t (E) (t 
on the south pole is replaced by a pseudospin-pseudospin 
structure Si{0) (g) S2{0) on T,{e). Note that the combi- 
nation only occurs for potentials with the same parity. 

Although the surface states solutions are stable in the 
presence of localized surface potentials that preserve T 
symmetry, two terms play an essential role in determining 
the energy position of Dirac point 



£bp—'H 



DP 



4m2(Ao 
4to2(Ao -I- 



Ae)K-A2 + A|) 



Ae)2 + (m2 



A§ 



A|) 



,(10) 



which implies that / and t^ potentials (Table IT| are able 
to tune the Dirac point from the midgap to the band 
edges ±m independently. 

Discussions. — The interactions of the topologically 
protected bands with surface potentials provides a ro- 
bust route to engineer and manipulate the topological 
surface states. In particular, an external surface poten- 
tial can raise or lower the Dirac point to the middle of the 
bulk gap and providing experimental access to the topo- 
logically protected band. This goal could be achieved 
by surface oxidation[28], or by other possible chemical 
processes |29f[3T] and interactions [32] on the surface. We 
also point out that the electrostatic gating[25, 5T1 1551 - 
155] . Ecs{r3)s{ can act to Stark shift the TI surface state 
into the gap, with field strength determined by the pen- 
etration length of surface states. Typically, the screened 
field ^ 100 meV-nm^^ are required for the surface states 
evanescent in a couple of quintuple layers. Our present 
results provide a framework to study the self-consistent 
band bending physics of real TI materials [HI [101 [^5H?7] 
and mean-field models of surface state many-body inter- 
actions. 

Our model is quite different from the fixed bound- 
ary condition (FBC) which arbitrarily clamps the sur- 
face state wavefunction to zero at the boundary. The 
FBC solution is insensitive to the mass inversion at the 
surface which topologically protects the surface bands, 
and consequently it provides no information about the 
energy of the symmetry protected degeneracy relative to 
the bulk bands or their interaction with surface-localized 
potentials. Furthermore FBC admits an infinite num- 
ber of (physically spurious) solutions with nonzero ener- 
gies for fci = fc2 = that satisfy an (incorrect) surface 
boundary condition. By contrast TBC guarantees that 
there is only one isolated solution, i.e., the Dirac point of 
surface bands protected by the change of bulk topology. 
Moreover, TBC demonstrates that the energy position of 



Dirac point is tunable in the bulk gap via the symmetry 
allowed scalar terms and surface potentials. 

On a non cleavage surface, dangling bonds and their 
reconstruction may add complexity to the surface state 
spectrum which can be modeled as surface potentials en- 
coded in the parameters Aq and Ae. The fingerprint 
of the novel spin texture near the Dirac point on non 
cleavage surfaces are determined by the bulk symmetries 
along with the topological stability of the surface spec- 
trum and may be accessible in ARPES and STM ex- 
periments. On an equatorial face, Zeeman coupling or 
magnetic disorder coupled to the spin degree of freedom 
do not generally open a gap, and the diamagnetic suscep- 
tibility is anticipated to be unusually anisotropic. More 
spectacularly, there is intrinsic charge redistribution in 
the surface bands near the corners of a TI that connect 
different crystal faces with intrinsically different Dirac 
point energies. 
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